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ABSTRACT 

Intensity  scintillation  in  a  laser  beam  at  0.63 
micrometers  in  the  marine  boundary  layer  has  been  studied 
over  a  4.3  kilometer  horizontal  path  across  Monterey  Bay 
and  also  from  shore  to  ship  at  San  Nicolas  Island.   Optical 
and  micrometeorological  measurements  of  the  refractive 

index  structure  constant,  CN,  agree  to  approximately  one 

_  g 
standard  devia^io^  giving  values  in  the  range  2.0  x  10 

-1  '3 
to  5.0  x  10   (meters)   '  .   Optical  and  meteorological 

data  correlation  improved  as  mean  wind  speed  increased. 

Evidence  was  found  to  support  the  reported  existence  of  a 

saturation  region  for  the  laser  beam  logarithmic  amplitude 

.variance . 
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I.   INTRODUCTION 

Atmospheric  effects  on  laser  beam  propagation  are 
divided  into  three  general  categories:   refraction,  absorp- 
tion, and  particulate  scattering.   This  investigation  was 
exclusively  concerned  with  atmospheric  refraction  effects. 
Refraction  may  be  further  subdivided  into  large-scale 
refraction,  where  the  beam's  spatial  path  is  altered,  and 
small-scale  ref^ac  ::ion,  where  the  beam  is  differentially 
focused  and  defocused.   Large-scale  refraction,  commonly 
referred  to  as  beam  wander,  was  found  to  be  negligible.   All 
research  efforts  were  directed  toward  the  small-scale  refrac- 
tion phenomenon  commonly  referred  to  as  beam  or  atmospheric 
scintillation. 

The (major  purpose  of  the  research  was  to  make  and  to 
correlate  direct  optical  measurements  of  laser  beam  scintil- 
lation with  meteorological  measurements  of  certain  atmospheric 
parameters.   Using  Rytov's  method  as  described  by  Tatarski 
[1]  ,  the  optical  and  meteorological  measurements  of  beam 
scintillation  were  compared  by  calculating  the  atmospheric 
refractive  index  structure  constant  which  statistically 
describes  atmospheric  refractive  index  fluctuations.   The 
research  team  attempted  to  predict  laser  beam  scintillation 
by  making  micrometeorological  measurements  of  atmospheric 
parameters . 

The  project  effort  consisted  of  two  primary  areas,  an 
optical  determination  of  the  refractive  index  structure 
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constant  (Cx,n)  an<3  a  meteorological  determination  of  the 
refractive  index  structure  constant  (CxjT) .   These  measure- 
ments were  made  simultaneously  and  compared  for  correlation. 
This  experiment  differs  from  previous  experiments  conducted 
by  Livingston  [2],  Ochs  [3]  and  Fitzmaurice  [4]  in  that  the 
laser  path  was  over  the  ocean  instead  of  over  land. 

This  research  project  was  funded  by  the  Naval  Ordnance 
Laboratory  and  was  divided  into  two  related  areas,  meteorol- 
ogy and  optics.   The  meteorological  research  group,  headed 

■ 

by  Dr.  Kennel*  Da/idson  of  the  Naval  Postgraduate  School, 
made  and  evaluated  the  micrometeorological  measurements  of 
atmospheric  turbulence.   The  optical  research  group, 
headed  by  Dr.  Eugene  Crittendon  (until  his  departure  for 
Australia  in  January)  and  later  by  Dr.  Alfred  Cooper,  made 
and  evaluated  the  optical  measurements  of  laser  beam 
scintillation.  _  i" , 

Due  to  the  wide  extent  of  each  of  the  areas  of  investiga- 
tion, each  group  worked  independently  of  the  other  to  obtain 
an  estimate  of.  atmospheric  turbulence.   These  estimates  were 
compared  for  correlation  after  data  reduction  and  analysis. 
The  inherent  nature  of  this  organization  promoted  speciali- 
zation of  the  personnel  involved  in  either  the  optical  or 
meteorological  areas  of  research.   The  author,  by  virtue  of 
his  background  in  electro-optics,  concentrated  his  efforts 
in  optical  research.   Hence,  this  report  is  devoted  pri- 
marily to  the  procedures  and  conclusions  of  the  optical 
research  group.   This  thesis  represents  a  comprehensive 
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description  of  the  efforts  and  results  of  the  optical 

propagation  study.  Meteorological  techniques  and  experi 

mental  results  are  described  briefly  when  they  relate  to 
optical  data. 


ir  i 
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II.   PROJECT  BACKGROUND 

A.   HISTORICAL  SUMMARY 

The  earliest  observations  of  scintillation  in  a  light 
beam  propagating  through  the  atmosphere  were  made  by  ancient 
astronomers  who  noticed  the  "twinkling"  of  stellar  images. 
Stellar  scintillation  was  measured  quantitatively  by  Middleton 
in  1949  [5]  and  1952  [6].   Chandrasekhar ,  in  1952  [7], 
estimated  the  phase  and  amplitude  fluctuations  of  starlight 
caused  by  atr jspheric  turbulence  by  means  of  a  geometrical 
optics  approach.   In  1953,  Mintzer  [8]  applied  the  Born 
approximation  to  the  scalar  wave  equation  to  describe  the 
propagation  of  sound  waves  through  a  turbulent  ocean,  a 
problem  similar  to  the  propagation  of  light  through  a  turbu- 
lent atmosphere. 

Muchmore  and  Wheelon  [9],  in  1955,  used  Chandrasekhar " s 
approach  to  study  line  of  sight  propagation  of  radio  waves 
in  the  troposphere.   They  calculated  the  variance  and 
spatial  correlation  functions  for  the  amplitude  and  the 
phase  of  the  propagating  wave.   Two  articles  by  Wheelon 
published  in  1959  [10,  11]  review  the  state  of  the  art  of 
wave  propagation  at  that  time.   These  articles  cover 
turbulence  theory' and  atmospheric  propagation  of  radio 
waves  which  is  a  similar  problem  to  propagation  of  optical 
beams. 

In  1960  and  1961,  the  Russian  monographs  by  Chernov  [12] 
and  Tatarski  [1]  were  translated  into  English.   These  works 
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described  the  method  proposed  by  Rytov  [14]  (see  Section 
III.D  for  the  derivation)  for  analyzing  the  diffraction  of 
light  by  ultrasonic  waves  and  later  applied  by  Obukov  [15] 
to  the  propagation  of  light  waves  through  a  turbulent 
atmosphere.   Strohbehn  [16]  presents  a  good  discussion  on 
the  range  of  validity  and  approximations  made  by  Rytov' s 
method.   The  calculations  in  this  project  were  based  on 
Tatarski's  formulation  of  Rytov's  method.   Rytov' s  method 
takes  into  consideration  multiple  scattering  and  diffrac- 
tion effects  .nd  thereby  represents  a  significant  advance- 
ment in  the  theory  of  atmospheric  propagation  of  planar  or 
spherical  electromagnetic  waves.   Strohbehn' s  article, 
published  in  1968,  is  an  excellent  review  of  atmospheric 
propagation  theory. 

Lumley  and  Panofsky  [17]  wrote  (in  1964)  a  monograph 
covering  many  facets  of  turbulence  theory  as  applied  spe- 
cifically to  the  atmosphere.   This  defined  the  -meteorological 
theory  used  by  the  Naval  Postgraduate  School's  research  group 
to  obtain  estimates  of  the  atmospheric  refractive  index  struc- 
ture constant.   In  1972,  Davidson,  in  a  Naval  Postgraduate 
School  report  [18],  comprehensively  described  the  meteoro- 
logical measurements  made  and  the  theory  applied  in  this 
project.   Articles  by  Wyngaard  (e_t  al . )  in  1971  and  Panofsky  [19, 

20]  in  1968  provide  further  references  on  meteorological 
turbulence  theory  as  applied  by  the  research  group. 
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B.   RECENT  DEVELOPMENTS 

In  his  monograph,  Wave  Propagation  in  a  Turbulent 
Medium  [1] ,  Tatarski  reported  the  results  of  some  of  the 
earliest  experiments  that  compared  micrometeorological 
estimates  to  optical  estimates  of  atmospheric  turbulence 
using  Rytov's  method.   These  experiments  were  performed  in 
the  late  1950s   using  an  incoherent  light  source.   The 
propagation  path  ranged  from  250  to  2000  meters  over  a 
Russian  steppe.   Experimental  results  strongly  supported 
his  formulation  of  propagation  theory.   The  fluctuations  of 
beam  intensity  and  phase  caused  by  atmospheric  turbulence 
were  found  to  have  a  log-normal  probability  distribution. 
Tatarski' s  description  of  the  scintillation  frequency  spectrum 
and  the  dependence  of  the  fluctuation  variance  on  path  length 
were  found  to  match  propagation  measurements. 

In  1966,  Schmeltzer  [22]  published  a  derivation  of  the 
covariance  of  the  logarithmic  intensity  and  phase  fluctua-  . 
tions  of  a  propagating  beam  of  light.   His  results  were 
more  general  than  Tatarski's  as  they  are  valid  for  colli- 
mated  and  focused  beams  as  well  as  plane  and  spherical 
waves.   However,  his  results  were  in  the  form  of  compli- 
cated double  integral  equations  and  so  were  not  directly 
useable.   Fried  evaluated  these  integrals  to  obtain  analytic 
expressions  for  the  logarithmic  intensity  variance  of  a  spheri- 
cal wave  [23]  and  for  a  collimated  or  focused  beam  [24],   His 
expression  for  spherical  wave  scintillation  was  in  agreement 
with  Tatarski's  results. 
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In  1969,  Ochs  [3]  repeated  Tatarski's  experiments  in 
Colorado  using  a  helium-neon  laser  as  a  light  source.   His 
experiment  found  several  deviations  from  Tatarski's  theory- 
including  the  existence  of  a  saturation  region  for  the  opti- 
cal measurements.   When  atmospheric  turbulence  was  very 
high,  Ochs  found  that  meteorological  measurements  would 
continue  to  measure  turbulence  values  but  the  optical 
measurements  would  saturate  at  a  lower  turbulence  value. 
Optical  and  meteorological  measurements  would  not  corre- 
late until  tl J  turbulence  had  decreased  below  the  satura- 
tion level  of  the  optical  measurements. 

Similar  experiments  were  performed  over  land  by 
Fitzmaurice  in  1970  [4]  and  by  Dowling  and  Livingston 
in  1971  [2] .   Both  groups  measured  the  scintillation  of 
0.6328  and  10.6  micrometer  laser  beams  and  compared  optical 
and  micrometeorological  estimates  of  atmospheric  turbulence. 
Both  groups  found  evidence  of  a  saturation  region. 

The  work  described  in  this  thesis  was  part  of  an  inter- 
departmental program  formed  as  a  combination  of  the  efforts 
of  existing  optical  and  meteorological  research  groups  at 
the  Naval  Postgraduate  School.   This  study  applied  Tatarski's 
formulation  of  Rytov's  method  to  describe  atmospheric  turbu- 
lence.  Optical  and  micrometeorological  measurements  of 
laser  beam  scintillation  and  atmospheric  temperature  fluctu- 
ations have  been  made  independently  and  then  analyzed.   Both 
types  of  measurement  yielded  values  for  the  index  of  refraction 
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structure  constant,  CN.   Theory  predicts  that  if  the 
turbulence  is  isotropic  and  homogeneous,  these  values  of  CN 
should  be  equal.   However,  two  major  unresolved  questions 
remained.   First,  does  saturation  of  optical  values  of 
CN  exist  in  the  high  turbulence  regime?   If  so,  what  param- 
eters determine  this  saturation?   Ochs  [3]  reported  that  the 
saturation  level  seemed  to  vary  with  mean  wind  speed  and 
the  time  of  day. 

Second,  Friehe,  of  Scripps  Institute,  reported  recently 
that  humidity  fluctuations,  which  were  neglected  by  Tatarski, 
may  contribute  significantly  to  refractive  index  fluctuations 
in  the  atmosphere  [25].   His  findings  have  not  been  general- 
ly accepted  by  propagation  researchers.   The  frequency  response 
of  available  humidiometers  is  inadequate  to  check  these  pre- 
dictions . 

Although  considerable  evidence  exists  to  support  Tatarski 's 
propagation  theory  over  land,  the  problem  of  propagation  in 
the  marine  environment  has  not  been  resolved.   If  humidity 
variations  influence  refractive  index  fluctuations,  laser 
propagation  over  the  ocean  may  not  be  described  by  Tatarski' s 
formulation  of  Rytov's  method.   This  thesis  examines  the 
validity  of  Tatarski' s  propagation  theory  in  the  marine 
boundary  layer. 
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III.   THEORY 

A.   GENERAL  REMARKS 

Atmospheric  turbulence  is  a  three-dimensional,  nonlinear, 
stochastic  process  that  occurs  in  air  whenever  certain  fluid- 
mechanical  parameters  are  exceeded.   Random,  dissipative 
fluctuations  in  the  homogeneity  of  the  atmosphere  will  cause 
the  velocity  profiles  of  the  air  to  vary  in  time  and  space. 
These  fluctuations  have  significant  effects  on  the  propaga- 
tion of  a  laser  beam  by  continuously  changing  atmospheric 
density.    This  causes  the  refractive  index  to  fluctuate 
with  time  and  space. 

Atmospheric  turbulence  (linear)  effects  on  laser  beam 
propagation  may  be  classified  into  two  categories.   Large- 
scale  inhomogeneities  will  uniformly  refract  the  laser  beam 
and  cause  the  beam  to  wander  in  space.   Small-scale  inhomo- 
geneities will  differentially  refract  the  laser  beam  and 
cause  interference  phenomena  or  scintillation.   Small-scale 
turbulence  alters  the  optical  path  length  over  different 
segments  of  the  beam  and  generates  interference  patterns 
between  adjacent  wave  fronts.   This  alters  the  initial  phase 
and  amplitude  distributions  of  the  propagating  beam  and 
destroys  spatial  and  temporal  coherence. 

Due  to  the  inherent  causes  of  turbulence,  it  is  impossible 
to  predict  the  transient  velocity  fields  that  arise  for  a 
given  set  of  boundary  flow  conditions.   Statistical  analysis 


20 


must  be  used  to  describe  atmospheric  turbulence.   By- 
analyzing  the  variations  in  atmospheric  homogeneity,  one 
may  describe  the  fluctuations  of  the  index  of  refraction 
over  a  given  path. 

Relating  the  fluctuations  in  a  propagating  laser  beam 
to  atmospheric  turbulence  is  a  three-part  problem.   The 
first  part  is  to  find  a  convenient  formulation  for  charac- 
terizing atmospheric  turbulence.   This  may  be  done  by 
using  statistical  estimators  of  the  covariance  and  the 
structure  function  to  describe  fluctuations  in  the  index 
of  refraction.   The  second  part  of  the  problem  is  to  derive 
the  equations  that  govern  the  propagation  of  electromagnetic 
waves  through  a  random  medium,  starting  with  Maxwell's 
equations.   These  equations  will  relate  the  magnitude  and 
phase  of  the  radiation  to  the  index  of  refraction  of  the 
optical  path.   The  final  part  of  the  problem  is  to  combine 
the  statistical  description  of  the  index  of  refraction  with 
the  propagation  equation.   The  result  will  be  an  analytic 
relationship  between  the  statistical  description  of  the 
refractive  index  and  a  statistical  description  of  the 
phase  and  amplitude  of  the  propagating  beam. 

The  entire  derivation  is  extremely  long  and  difficult. 
Consequently  only-  the  physically  most  important  points  of 
the  derivation  are  discussed  here,  and  much  of  the  turbu- 
lence theory  and  mathematics  are  omitted.  A  complete 

treatment  of  the  propagation  problem  may  be  found  in 
references  1,  17,  22,  23  and  24. 
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In  order  to  minimize  the  effects  of  beam  wander  on  the 
scintillation  measurements,  the  laser  beam  was  diverged. 
The  propagating  diverged  beam  was  approximated  by  a  spheri- 
cal wave.   Thus  the  general  propagation  equations  obtained 
were  evaluated  for  the  case  of  a  spherical  wave. 

B.   TURBULENCE  THEORY 

Wave  propagation  in  a  random  continuous  medium  is 
governed  by  a  stochastic  wave  equation,  whose  index  of 
refraction,  n(r,t),  is  a  continuous  random  point  function, 
and  characterizes  the  transmission  medium.   For  line  of 
sight  propagation  on  the  earth,  optical  paths  are  relatively 
short  enabling  one  to  say  that  the  transit  time  of  a  given 
phase  front  between  the  source  and  receiver  is  much  less 
than  the  time  constant  associated  with  changing  spatial 
distribution  of  the  refractive  index  (typically  10   seconds) 
The  time  dependence  of  n(r,t)  will  be  accounted  for  by 
letting  each  n(r,t.)»  for  i  =  1,  2,  3...  be  considered  a 
single  realization  of  the  random  process,  n(r).   Hence,  one 
can  characterize  the  transmission  medium  by  a  refractive 
index  in  the  form 

CD         n(f)  =  <n(r)>  +  n^r)  =  1  +  n^r) 

where  n,  (r)  represents  the  fluctuations  in  n(r), 

Tatarski  [1]  assumed  that  the  power  spectrum  of  the 
turbulence,  §_ (K) ,  can  be  described  by  the  empirical 
Kolomogorov  model  where 
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(2)       *n(K)  =  0.033  CN2  K'11/3  exp(-K2/Km2) 


C»T  is  a  constant  that  parameterizes  the  total  amount  of 
N 

energy  in  the  turbulence.   Tatarski  showed  that  CN  can  be 
related  to  a  statistical  estimator  by 


(3)  Dn(d)  =  CN2  d2/3 


where  D  (d)  is  the  structure  function  of  the  refractive 
n 

index  given  by 


Dn(d)  -  <[n(r)  -  n(r+d)]2> 

and  hence 

m              r  2  -  <[n(r)  -  n(r+d)]2) 
•C4J  CN ^273 


Assuming  isotropic,  homogeneous  turbulence,  the  functional 

arguments  r  and  d  can  be  changed  from  vectors  to  scalar 

distances. 

The  quantity  K,  the  spatial  wave  number,  is  given  by 

K  =  2tt/£  where  I   is  commonly  interpreted,  in  some  sense, 

as  the  size  of  the  turbulent  eddy.   This  turbulence  spectrum 

is  valid  for  values  of  K  that  are  in  the  inertial  subrange, 

that  is,  the  range  where  £Q  <  £  <  Ln.   The  inner  scale,  lQ, 

characterizes  the  size  of  turbulent  eddies  where  viscous 

dissipation  of  turbulent  energy  becomes  predominant  (for 

values  of  1   less  than  £A) .   K   is  the  spatial  wave  number 

0     m         r 

corresponding  to  the  inner  scale  (i.e.,  K  =  2-n/iQ)  .   Values 
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of  I   greater  than  the  outer  scale,  Ln,  characterize  the 
region  of  the  turbulence  spectrum  where  the  eddies  are  not 
homogeneous  and  where  energy  is  added  to  the  turbulence 
spectrum.   For  values  of  I   in  the  inertial  subrange,  the 
motion  of  the  larger  turbulent  eddies  causes  the  formation 
and  transfer  of  energy  to  smaller  eddies  with  very  little 
viscous  dissipation.   These  eddies  then  form  and  trans- 
fer energy  to  smaller  eddies  and  so  on  until  eventually, 
at  approximately  I   =  &n,  viscous  effects  become  the  pre- 
dominant method  of  energy  dissipation  and  turbulent  energy 
is  converted  to  thermal  energy.   Typically,  for  atmospheric 
turbulence  in  the  marine  boundary  layer,  l„    is  of  the  order 
of  a  few  millimeters  and  L0  is  of  the  order  of  a  few  meters 
(depending  on  the  stability  conditions  and  height  above 
the  boundary) . 

Turbulence  theory  shows  that  the  covariance  of  the 
refractive  index  for  homogeneous,  stationary  conditions 

(5)  '  Cn(d)  =  <n(r+d)  n(r)> 

is  the  three-dimensional  Fourier  transform  of  the  turbulence 
energy  spectrum 

(6)  C  (d)-  =  47r/d  /°°*  (K)  K  sin(Kd)  dK  . 

n  o   n 

Also,  the  structure  function  of  the  index  of  refraction, 
given  by  equations  (3)  and  (4) ,  is  related  to  the  turbulence 
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energy  spectrum  by  the  transform  relation 

(7)  Dn(d)  =  8tt  /°°  a    -    SiKdKd)^n(K)  K2  dK 

under  conditions  of  isotropy,  homogeneity  and  stationarity . 
Given  equations  (2)  and  (6) ,  one  can  relate  C  (r)  to  C^. 
This  means  that  C^  can  be  used  to  describe  the  statistics 
of  the  fluctuations  of  the  refractive  index  in  the  atmosphere 

C.   DERIVATION  OF  THE  SCALAR  WAVE  EQUATION 

The  equations  describing  the  propagation  of  an  electro- 
magnetic wave  in  the  atmosphere  are  derived  from  Maxwell's 
equations  for  a  non-conducting  medium  without  charge  sources. 

(8)  (a)    V  •  D  =  0 


00 

V    •    B    =    0 

(c) 

»xfi-  -ff 

(d) 

V   XH   =|| 

with  the  additional  relations 

(e)  D  =  eE       with  e  =  e(r) 

(f)  B  =  yH       with  y  =  y   =  constant. 
Taking  the  curl  of  equation  (8c)  one  gets 

(9)      V  X  V  X  E  =  V(V  •  El  -  V2  E  =  yV  X  |§   . 

a  t 

Assuming  a  time  dependence  of  e"la)   for  all  field  quantities 
and  substituting  equation  (8d)  into  the  right  hand  side  of 


25 


equation  (9) ,  one  finds 


do)   v(v  •  I,  -  v^e  -  -„  KpLHi  ,  _v  |_  (|flj 


2- 

ye  — 9 
8tZ 


2  £ 

=  U£0)   E  . 


Starting  with  Maxwell's  equation  (8a)  and  using  a  well  known 
identify  for  the  divergence  operator, 

(11)     V«D=V»eE=eV»E+E»VG=0 

or 

{12)  V  •  E  =  -  £  'JZ      . 

Substituting  equation  (12)  into  equation  (10)  gives 

(13)  „2  e      2  6    ^rE  •  Ve. 

K      J  V   E  +  yew   E  =  -V(  )  . 

The  index  of  refraction  for  a  medium  is  given  by 

(14)  n  =  £  =  c(ye)1/2  =  (const) (e) 1/2 

where  c  equals  the  velocity  of  propagation  in  free  space, 
V  equals  the  velocity  of  propagation  in  the  medium  and  the 
wave  number  is  given  by 

(15)  k  =  |-\ 
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Hence 


2 
(16)  a)  ye  =  (^y)  (c  ye)  =  k   n      and 

c 


(17)  V2  E  +  k2  n2  E  =  -V(E  '  ^E)  . 


2 
Realizing  from  equation  (14)  that  n   =  (const) e,  then 

-  2 
—  =  -^—   =  2Vln(n)  and  equation  (17)  becomes 
e  n 


(18)         V2  E  +  k2  n2  E  =  -2V[E  •  Vln(n)]  . 


The  term  -2V[E  •  Vln(n)]  couples  the  various  components 
of  the  field  together  and  is  responsible  for  polarization 
fluctuations  in  the  linearly  polarized  beam.   Fitzmaurice  [4] 
shows  that  if  the  deviations  of  the  refractive  index  from 
its  mean  value  is  very  small  (typically  n-./(n>  ~  10   )  and 
if  the  wavelength  of  the  propagating  wave  is  much  less  than 
the  inner  scale  (&q)  of  the  turbulence,  then  one  can  assume 
that  the  polarization  fluctuation  term  is  negligible.   This 
gives  the  scalar  wave  equation 


(19)  V2  E  +  k2  n2  E  =  0 


where  E  is  any  rectangular  component  of  E(r). 
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D.   RYTOV'S  METHOD 

At  this  point,  Rytov  [14]  transformed  equation  (19)  to 


(20)         V2ln(E)  +  |Vln(E)|2  +  k2  n2  =  0  . 


Now  let 

(21)  E(r)  =  A(r)  exp[iS(r)] 

and  from  equation  (1),  n(r)  =  1  +  n, (r) . 
Define  the  complex  phase,  ip,  as 

(22)  E(r)  =  exp  |>(r)  ]  =  exp  [lnA(f)  +  iS(r)] 
and  therefore, 

(23)  ln(E)  =  ln(A)  +  iS  =  ^(r). 
Substituting  this  result  into  equation  (20)  gives 


(24)  V2ip  +  |V^|2  +  k2(l  +  nx)2  =  0 


Now  set 


(25)  ip  =  ip   +  \p1        so  that 


(26)  E  =  exp(4>Q  +  ip1)  =  Eq  exp^)  . 

E  is  the  unscattered  wave  which  satisfies  the  equation 

(27)  V2E   +  k2  E   =  0  J 

oo 
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or  the  transformed  equation 


(28)         V2ln(E  )  +  |Vln(Eo)|2  +  k2  =  0 


Substituting  equation  (22)  into  equation  (28) , 


(29)  y2\p      +  |Vip  I2  +  k2  =  0  . 

Now,  substituting  equation  (25)  into  equation  (24)  gives 


(30)      V2^      +    V2ip1    +    |V^q    +    V^1|2    +   k2(l    +    2nx    +   n^)    =    0 


and   subtracting    this    from   equation    (29)    gives    after   some 
manipulation 

(31)      V2ipx    +    |V^1|2    +    2V^Q    •    V^1    +   k2(2nx    +   i^2)    =    0. 


This    is    the   nonlinear   Ricatti    equation.       Since   n,    is   much 

2      2  2 

less    than  one,    k  n,      is   much   less    than   2k  n-,    and  may  be 

neglected.      Similarly,    since    |Vif>.  |    <<    |  Vi/>    |    ~    k,    then 

2 
|  Vip -.  |      is   much   less,  than   ZVifi      •    V^      and   one   can  also 

i  2 
neglect   the    |  V\p  ,|      term.      Finally,    one   can  write   equation 

(31)  as 

(32)  72^1    +    2V\f;o    •    Vi^    +    2k2    i^    =    0    . 


The  assumption  that  |v>- |  <<  \Vty    |  .  k  is  the  only 
significant  assumption  made  in  Rytov's  method.   It  is 
equivalent  to  saying  that  the  change  in  the  variation 
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of  the  complex  phase,  iK  ,  over  a  distance  of  about  one 
wavelength  in  much  less  than  2tt  radians.   Realizing  that 
the  wavelength  is  much  less  than  the  inner  scale,  £„,  it 
is  difficult  to  imagine  the  phase  changing  significantly 
in  a  propagation  distance  less  than  l^.      This  approximation 
is  generally  referred  to  as  the  Rytov  approximation  in  the 
literature. 

A  new  quantity  is  now  defined 


r 


(33)  W  =  ^e  °  =  i|;1Eq 


so  that 


(343  *i  -  H-  -  "  e"*° 

o 


Substituting  equation  (34)  into  equation  (32)  and  expanding 
terms  gives 

(35)  V2W  +  k2W  =  -2k2n1E   . 

1  o 

This  differential  equation  has  the  solution 


(36)    W(r)  =  /   2k2n,(r')  E  (r«)  G(f  -  r')  dV» 


where  G(r  -  r')  is  the  Green's  function 


(xn\  rrZ.        ;n         exp[ik(r    -    r')]         £  .     . 

C37J  G(r   -    r')    =   — r  L    _  "- — J-±-       for   a  point   source. 

4tt  (r   -    r') 
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Therefore, 


(38)   iK(f)  =  |~  =  i-f-.  7  2  k2  n  (r')  EQ(r')  G(f-  f')  dV 
o    oK    J    V 


Equation  (38)  with  n, (r)  as  a  random  continuous  function 
characterizes  the  scintillation  in  the  propagating  wave  if 
n-,(r)  is  known  analytically. 

It  is  not  practical  to  measure  the  index  of  refraction 
at  every  point  over  the  optical  propagation  path.   Instead 
it  is  necessary  to  statistically  characterize  n, (r)  using 
the  covariance  function,  C  (r) .   This  brings  one  to  the 
third  part  of  the  problem.   The  statistics  of  \p-.  (r)  must 
be  calculated  starting  from  equation  (38)  and  the  statisti- 
cal description  of  the  covariance  of  the  refractive  index 
as  given  by  equation  (6) . 


E.   APPLICATION  OF  TURBULENCE  THEORY  TO  THE  PROPAGATION 
EQUATION 

After  making  many  transformations,  substitutions  and 

integrations,  Schmeltzer  [22]  obtained  the  statistical 

2 
expression  for  the  logarithmic  amplitude  covariance,  C-.(d) 

for  a  beam  propagating  in  a  random  medium  which  Fried  [23] 

reduced  to 


^  C£A^  -  {£  ^  {  /  ds  /\(K1/2)  Jo(K1/2d  f) 


X  [1  -  exp  (KS^2"  S))]  dK  } 
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for  the  case  of  a  spherical  wave.  The  quantity  Mz"  is 
the  optical  propagation  path  length  and  "s"  is  a  dummy 
variable  of  integration  with  the  dimensions  of  length. 

Fried  integrates  this  expression  for  a  spherical  wave 
over  s  and  K  to  obtain  the  relation 


(40)    C^A(d=0)  =  C£2  =  0.124  k7/6  z11/6  CN2   where 


(41)  C£A2  =  ([In (A)  -  ln<A)]2> 

and  A  is  the  amplitude  of  the  wave  function,  |E(r)|. 

The  logarithmic  intensity  variance,  Cn-r2,  is  a  more  readily 
measurable  quantity 

(42)  Cn2    =     ([ln(I)    -    ln<I}]2> 

2 
where  I  is  the  intensity  of  the  radiation  (I  =  A  ).   Thus, 

2    /ro,  „.<»    „,   /,x  -,2> 


(43)  C£IZ  =  ([21n(A)  -  21n  <$  ]£) 


or  finally 


=  4  ([ln(A)  -  ln<&]2) 


=  4  C£A2   and 


(44)        •  C£T2  =  4[0.124  k7/6  z11/6  CN2] 


(45)         CN  =  1.42  k"7/12  z"11/2  C£I 


The  foregoing  is  a  brief  summary  of  the  major  steps  in 
the  development  of  the  governing  equation  for  the  propagation 
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statistics  of  a  spherical  wave.   Much  of  the  detailed 
mathematics  involved  in  the  derivation  has  been  omitted. 
For  the  complete  derivation,  the  original  papers  by 
Tatarski  [1] ,  Schmeltzer  [22],  and  Fried  [23]  should  be 
consulted.   However,  the  approximations  involved  in 
reducing  the  equations  to  analytic  form  are  interesting 
and  will  be  analyzed  here. 

F.   THEORETICAL  ASSUMPTIONS 

The  major  assumption  that  underlies  the  entire  theory 
is  that  the  turbulence  energy  spectrum  can  be  described  by 
the  empirical  Kolomogorov  model  [equation  (4)].   This  model 
was  developed  in  connection  with  wind  tunnel  experiments. 
Meteorological  measurements  by  Tatarski  [1] ,  by  Wyngaard, 
Izumi  and  Collins  [19]  and  by  Panofsky  [20]  verified  that 
it  was  applicable  in  the  atmospheric  boundary  layer  over 
the  inertial  subrange^  It  was  also  implicitly  assumed  that 
only  those  turbulent  eddies  in  the  inertial  subrange  contribute 
significantly  to  beam  scintillation.   Eddies  with  £  >  Lfi  would 
tend  to  cause  beam  wander.   Eddies  with  I   <   Lfl  have  little 
effect  on  the  scintillation  because  viscous  forces  cause 
rapid  collapse. 

The  assumption  that  A  <<  l~   was  used  to  justify  ignoring 
the  cross  polarization  term  in  the  scalar  wave  equation 
(equation  18) .   The  helium-neon  laser  used  in  the  experiment 
had  a  wavelength  of  0.6328  micrometers  which  is  much  less 
than  the  dimensions  of  the  inner  scale  (£n  was  of  the  order 
of  a  few  millimeters). 
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In  order  to  use  the  perturbation  technique  to  express 
the  variations  in  the  index  of  refraction,  the  fluctuations 
in  refractive  index  were  assumed  to  be  much  smaller  than  the 
mean  value.   Meteorological  measurements  confirm  that  the 
index  of  refraction  in  the  atmosphere  does  not  increase 
by  more  than  a  few  per  cent  due  to  density  variations; 
hence,  this  assumption  seems  justified. 

The  assumption  that  A  z/£f)  <<1  is  a  very  weak  restriction 
at  optical  wavelengths.   For  a  helium-neon  laser  and  an 

inner  scale  value  of  one  millimeter,  the  propagation  path 

8 

would  have  to  be  greater  than  approximately  10   meters  for 

this  assumption  to  become  invalid. 

The  final  assumption  was  that  turbulence  was  isotropic 
and  homogeneous.   This  approximation  is  valid  for  large 
areas  with  a  smooth  regular  topography.   Ocean  conditions 
are  conducive  to  isotropic  homogeneous  turbulence.   This 
approximation  means  that  a  point  turbulence  measurement  can 
characterize  the  turbulence  over  an  entire  propagation  path. 

The  experimental  conditions  seemed  to  satisfy  the  require- 
ments of  Tatarski's,  Schmeltzer's  and  Fried's  theory  for  the 
propagation  of  a  spherical  wave  in  a  random  medium.   A  goal 
of  this  experiment  was  to  test  the  validity  of  their  propa- 
tion  theory  in  the  marine  boundary  layer. 

G.   PROBLEM  APPLICATIONS 

The  calculation  of  C*.  from  equation  (45)  required  values 
of  the  wavelength  of  the  laser  beam,  the  optical  path  length 
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and  the  variance  of  the  logarithmic  intensity  fluctuations 
of  the  laser  beam.   The  wavelength  of  the  laser  was  known 
to  be  0.6328  micrometers.   The  optical  propagation  path 
length  was  measured  on  a  nautical  chart.   Log-intensity 
variance  was  measured  using  the  laser  scintillometer  which 
is  described  in  the  equipment  section  of  this  thesis.   The 
calculated  values  of  C„  were  designated  C„0,  indicating  that 
they  were  derived  from  optical  measurements.   These  values  of 
CN0  were  compared  with  values  of  CN  that  were  calculated  from 
meteorological  measurements  (designated  C„T) .   A  brief 
explanation  of  the  theory  and  equations  that  describe  the 
calculation  of  C„T  is  appropriate  at  this  point.   The 
meteorological  theory  is  described  in  detail  by  Panofsky 
•[20]  and  by  Wyngaard,  Izumi  and  Collins  [19]. 

H.   METEOROLOGICAL  TURBULENCE  MEASUREMENTS 

Variations  in  the  index  of  refraction  in  the  atmosphere 
are  caused  by  variations  in  the  density  of  the  air.   These 
variations  are  due  to  fluctuations  in  the  temperature, 
pressure,  humidity  and  perhaps  gravity  waves.   Empirical 
measurements  [4]  show  that  the  index  of  refraction  at  opti- 
cal wavelengths  can  be  represented  by  the  equation 

(46)     n(r)  =  1.0  +  77.6  ££|i  (1  +  °-°°75)  X  10'6 

A. 

where  atmospheric  pressure,  P(r)  is  in  millibars,  absolute 
temperature,  T(r),  is  in  Kelvins,  and  wavelength,  A,  is  in 
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micrometers.   Evaluation  of  equation  (46)  at  a  wavelength 
of  0.6328  micrometers  gives 

(47)  n(r)  =  1.0  +  79.0  ^[||  X  10'6  .   . 

It  is  generally  accepted  in  atmospheric  propagation 

i 
studies  that  temperature  fluctuations  are  the  predominant 

mechanism  causing  index  of  refraction  fluctuations  in  the 

1  2 

atmosphere.    Equation  (4)  defines  C^  as 

m             r  2  _  ([n(r)  -  n(r  +  d)]2) 
(4)  LN   -        -^3 . 

2 
One  can  also  define  a  temperature  structure  constant,  CT  by 

(48)  CT2  =  ([TW  -J(r+d)]2)  _ 


Rearranging  terms  gives 

([(79.0  - X  10"6)  -  (79.0  - X  10'6)]2\ 

(50)  c  2  „  ]_  T(r)  T(r+D)  ' 


'N 


d2/3 


and  assuming  that  the  temperature  variations  are  much  smaller 
than  the  mean  temperature,  then 

lbiJ     •     TIrT  "  T(r+d)  "  TJFY   ~  T     and 


Recently,  Friehe  of  Scripps  Institute,  in  a  private 
communication,  suggested  that  humidity  variations  may  contrib- 
ute significantly  to  refractive  index  fluctuations  under  con- 
ditions of  high  humidity.   This  has  not  been  generally  accepted, 
It  is  difficult  to  measure  humidity  fluctuations  due  to  a  lack 
of  humidiometers  that  will  measure  fluctuations  with  a  fre- 
quency greater  than  15  Hz  [25], 
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((79. 0  L-  X  10"6)2  [T(r)  -  T(r+d)]2\ 

2    '      T  ' 

(52)    CN  -^ 

a 
Substituting  equation  (48)  into  equation  (52)  gives 


(53)  CN2  =  (79.0  ~   X  10"6)2(CT2)  . 


Thus  the  analytic  expression  for  C„T  is 


(54)  CNT  =  (79.0  ^  X  10"6)  C, 


In  order  to  calculate  C„T  using  equation  (54)  measurements 
of  mean  pressure,  mean  temperature  and  CT  as  given  in  equation 
(48)  must  be  made.   These  measurements  will  be  described  in 
the  experiment  section  of  this  thesis. 


37 


IV.   EXPERIMENTAL  APPARATUS  AND  PROCEDURES 

This  section  will  describe  the  apparatus  and  procedures 
used  to  make  measurements  during  the  project.   The  specific 
results  and  conclusions  derived  from  the  measurements  will 
be  reported  in  later  sections. 

Two  major  systems  were  used  to  measure  the  data  obtained 
during  the  three  main  experiments.   The  scintillometer 
system  measured  the  standard  deviation  of  the  logarithmic 
intensity  distribution,  C?T,  for  the  transmitted  laser  beam. 
These  measurements  furnished  the  data  for  the  calculation  of 
CNn  from  equation  (45).   The  thermosonde  system  measured  the 
temperature  structure  constant,  CT,  mean  temperature,  T, 
and  the  mean  pressure,  P.   These  measurements  furnished 
the  data  used  in  the  calculation  of  CNT  from  equation  (54) . 
Additional  experiments  were  performed  to  measure  the  intensity 
stability  of  the  laser,  the  power  spectrum  of  the  scintilla- 
tion and  the  degree  of  truncation  of  the  laser  beam  by  the 
laser  output  optics.   These  experiments  verified  the 
scintillometer  system's  performance  and  accuracy.   The 
equipment  used  in  the  experiments  will  be  subsequently 
described  in  this  section. 

A.   THE  SCINTILLOMETER  SYSTEM 

The  scintillometer  system  measured  the  logarithmic 
intensity  variance  of  the  transmitted  laser  beam.  The 
main  components  of  the  system  are  illustrated  in  Figure  1. 
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A  low  powered  helium-neon  laser  transmitted  a  beam  to  an 
optical  receiver.   The  receiver  detected  the  incoming  beam 
and  generated  an  electrical  voltage  proportional  to  the 
instantaneous  intensity  of  the  beam.   This  signal  was  ana- 
lyzed by  the  laser  scintillometer  and  recorded  on  magnetic 
tape.   The  standard  deviation  of  the  logarithmic  intensity 
of  the  laser  beam  was  computed  by  analogue  processing  methods. 
This  section  analyzes  in  detail  the  operation  of  this  system. 

The  laser  output  and  pointing  system  controlled  the  spot 
size  and  direction  of  the  beam  (Figure  2) .   The  laser  beam 
was  focused  by  the  first  lens  to  a  position  near  the  focal 
point  of  the  output  lens.   An  adjustable  mount  was  pro- 
vided to  give  the  output  lens  two  degrees  of  translational 
freedom.   The  spot  size  was  controlled  by  changing  the  focal 
point  of  the  output  lens  with  respect  to  the  fixed  focal 
point  of  the  first  lens.   During  normal  operation,  the  beam 
was  diverged  to  minimize  beam  wander  effects.   A  telescope 
with  cross  hairs,  bore-sighted  with  the  laser,  provided  an 
accurate  way  of  aiming  the  beam.   Beam  alignment  with  the 
target  was  controlled  by  elevation  and  azimuthal  adjustment 
screws.   A  chopper  was  provided  to  interrupt  the  beam  for 
ten  seconds  of  every  minute  during  which  adjustments  were 
made  on  the  background  compensator  at  the  receiver. 

After  transmission  through  the  atmosphere,  a  portion  of 
the  scintillating  beam  plus  background  radiation  was  collected 
by  the  input  optics  in  the  detector  package  (Figure  3) .   The 
incoming  radiation  was  passed  through  a  0.63  micrometer 
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bandpass  filter  and  focused  by  a  lens  through  two  small 
apertures.   The  detector  package  was  aligned  in  such  a 
manner  as  to  pass  the  laser  beam  plus  background  radiation 
through  the  first  aperture  to  the  first  silicon  photocell. 
The  second  aperture  passed  only  the  background  signal  which 
was  reflected  from  a  small  prism  to  the  second  photocell. 
The  photocells  generated  an  output  voltage  proportional  to 
the  intensity  of  the  radiation  that  impinged  upon  the  photo- 
sensitive surfaces. 

The  two  signals  were  then  subtracted  in  a  differential 
amplifier  (Figure  4)  thus  cancelling  the  background  signal 
and  passing  only  a  signal  that  was  linearly  proportional  to 
the  laser  radiation.   During  operation,  the  laser  beam  was 
interrupted  for  ten  seconds  every  minute,  during  which  time, 
the  background  signals  from  the  two  photocells  could  be 
electronically  adjusted  in  order  to  exactly  cancel  each  other. 

The  scintillation  signal  (Figure  5)  was  amplified  by  a 
factor  of  five  and  then  passed  through  a  logarithmic  ampli- 
fier.  The  logarithmic  amplifier  computed  the  logarithm  of 
the  original  signal  and  sent  it  to  a  wide  band  power  ampli- 
fier.  The  signal  level  was  multiplied  by  a  factor  of  ten 
and  recorded  on  magnetic  tape  with  an  F.M.  conversion 
instrumentation  tape  recorder.   The  signal  was  then  fed  to 
a  multi-channel  pulse  height  analyzer  which  sampled  the  signal 
level  periodically.  The  pulse  amplitude  was  used  as  a 
criterion  to  address  one  set  of  200  analyzer  memory  channels. 
These  channels  counted  the  number  of  pulses  of  a  given 
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voltage  level  that  were  accumulated  over  a  variable  time 
period  (usually  around  four  minutes).   The  pulse  height 
analyzer  output  consisted  of  a  CRT  display  of  the  number  of 
pulses  of  a  given  signal  level  as  a  function  of  the  logarithm 
of  the  beam  intensity.   This  corresponded  directly  to  a  non- 
normalized,  probability  distribution  of  the  logarithm  of 
the  laser  beam  intensity,  ln(I).   The  CRT  display  was 
photographed  using  a  polaroid  camera  and  the  photographs 
(example  shown  in  Figure  6)  were  later  analyzed  to  determine 

CN0* 

To  aid  in  the  analysis,  a  calibration  signal  generator 

provided  direct  current  calibration  signals  at  25,  50, 
100  and  200  millivolt  levels  which  were  accumulated  and 
.stored  on  a  second  set  of  200  pulse  height  analyzer  channels. 
The  two  sets  of  200  channels  could  be  displayed  in  such 
a  manner  as  to  superimpose  the  calibration  signals  and  the 
scintillation  pulse  height  distribution.   Amplifier  and 
pulse  height  analyzer  memory  drift  were  found  to  be  negli- 
gible for  the  system. 

Measurements  by  Tatarski  [1],  Ochs  [3]  and  Dowling  [4] 
have  shown  that  the  probability  distribution  of  the  logarithm 
of  beam  intensity  is  Gaussian.   This  was  expected  as  a 
direct  result  of  the  central  limit  theorem  and  the  randomness 
of  atmospheric  turbulence.   Realizing  that  the  pulse  height 
distribution  was  Gaussian,  then  the  maximum  point  corresponded 
to  the  mean  value  of  the  distribution  and  the  standard  devia- 
tion of  the  logarithmic  intensity,  C„T,  could  be  measured 
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directly  from  the  photographs.   Realizing  that  for  a 
Gaussian  distribution,  the  full  width  at  half  maximum,  Xpi/2, 
was  related  to  the  standard  deviation  by 

(56)  Cn    =  0.425  Xpl/2 

then  the  standard  deviation  could  be  determined  from  the 
following  equation. 

(57)  C^    =  0.425  S(Xpl/2  in  cm) 

"S"  is  a  scaling  parameter  that  related  the  distance  along 
the  abcissa  in  centimeters  to  the  signal  value  in  log- 
intensity  units.   S  was  determined  by  measuring  the  distance 
in  centimeters  between  calibration  marks  and  by  applying 
the  relation 

[ln(V2)  -  lnO^)  .] 

(58)  S      (X9cm  -  X,cm) 


v2 

In  ^ 
vl 


(AX)  cm 

where  V  is  the  calibration  voltage,  X  is  the  distance  along 
the  abcissa  in  cejitimeters  and  AX  is  the  distance  in  centi- 
meters between  two  calibration  marks,  V,  and  V-.   Once  Cpj 
was  determined,  Cj.  was  computed  using  equation  (45) . 

Implicit  in  the  system  computation  of  C5T  was  the 
approximation  that  In  (l)    =  (ln(I)^>.   This  approximation 
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was  tested  by  computing  both  quantities  for  typical 
probability  distributions  of  intensity  and  comparing  the 
difference.   The  error  introduced  for  large  values  of  CN~ 
was  less  than  2%   and  was  correspondingly  smaller  for  lower 
values  of  C„n.   This  error  was  judged  negligible  compared 
to  the  experimental  error  involved  in  measuring  CNT  and 
was  ignored. 

The  dynamic  range  and  frequency  response  of  the 
scintillometer  were  wide  enough  to  process  all  scintilla- 
tion signals  without  distortion.   The  logarithmic  amplifier 
had  a  dynamic  range  of  80  db  and  distorted  voltage  levels 
from  the  differential  amplifier  that  were  less  than  one 
millivolt  or  greater  than  ten  volts.   By  choosing  an 
appropriate  gain  setting,  the  scintillation  signal's  mean 
value  at  the  differential  amplifier  output  terminal  was 
controlled  to  be  approximately  100  mV.   All  other  system 
components  had  dynamic  ranges  greater  than  80  db .   The 
minimum  signal  levels  detected  during  the  experiments  were 
around  five  millivolts  while  the  maximum  levels  were  below 
one  volt.   System  frequency  response  was  flat  from  dc  to 
1.5  KHz.   However,  the  scintillation  power  spectrum  was 
measured  to  be  negligible  beyond  300  Hz  and  hence,  the 
signal  frequency  response  was  undistorted. 

B.   THE  THERMOSONDE  SYSTEM 

The  meteorological  equipment  used  for  measuring  C„T 
consisted  of  a  platinum  wire  thermosonde  unit,  a  barometer, 
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a  thermometer  and  a  psychrometer .   The  platinum  wire 
thermosonde  unit  measured  the  temperature  structure  function, 
CT,  by  measuring  the  changes  in  the  difference  of  resistance 
between  two  exposed  platinum  wires  mounted  on  probes.   The 
resistance  of  each  wire  was  proportional  to  the  temperature 
of  the  wire.   The  system  electronics  used  analogue  processing 
techniques  to  compute  CT  from  the  equation 

r  2  _  ([T(r)  -  T(r+d)]2) 
^T  ^273 

where  the  value  of  "d' ,  the  probe  separation  distance,  was 
constant.   This  system  had  a  linear  frequency  response  from 

dc  to  1  KHz  and  a  wide  dynamic  range.   The  system  generated 

2 
a  voltage  signal  corresponding  to  CT  which  was  ted  to  a 

strip  chart  recorder.   The  thermosonde  output  was  recorded 

as  a  function  of  time. 

Mean  pressure  was  recorded  as  a  function  of  time  from 

a  barometer.   Mean  temperature  was  read  and  recorded  from  a 

mercury  thermometer.   Relative  humidity  was  also  read  and 

recorded  from  a  psychrometer .   CNT  was  calculated  from 

measurements  of  CT,  mean  temperature  and  mean  pressure 

using  equation  (54) . 

C.   SCINTILLOMETER  PERFORMANCE  VERIFICATION  EXPERIMENTS 

Three  minor  experiments  were  performed  in  conjunction 
with  the  scintillation  measurements  to  determine  the 
scintillation  power  spectrum,  laser  output  intensity 
stability  and  beam  truncation  by  the  output  optics. 
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Each  experiment  was  necessary  to  verify  the  performance  of 
the  scintillometer  system. 

1.   Frequency  Response  and  Dynamic  Range 

In  order  to  assure  that  all  laser  beam  signal  levels 
were  being  accumulated  by  the  pulse  height  analyzer,  the 
dynamic  range  and  frequency  response  of  the  system  were 
tested.   The  system  response  to  a  wide  range  of  input  signal 
levels  was  found  to  be  logarithmic  for  signals  greater  than 
one  millivolt  and  less  than  ten  volts.   Typical  log-intensity 
distributions  ranged  from  approximately  15  to  400  millivolts 
and  all  measured  distributions  were  within  the  system's 
dynamic  range. 

Similarly,  frequency  response  tests  showed  that  the 
system  would  record  and  process  signals  with  a  frequency 
range  of  dc  to  1.5  KHz  without  distortion.   It  was  then 
necessary  to  determine  the  scintillation  power  spectrum  to 
justify  the  frequency  response  of  the  scintillometer  system. 
An  analogue  Fourier  spectrum  analyzer  was  designed  and 
assembled  using  a  wave  analyzer,  calibration  oscillator  and 
an  X-Y  plotter  (Figure  7) . 

The  FM  instrumentation  tape  recorder  fed  a  measured 
scintillation  signal  through  a  power  amplifier  to  a  wave 
analyzer.   The  wave  analyzer  passed  the  signal  through  a 
tunable  bandpass  filter  (bandwidth  of  3  Hz)  and  then 
generated  a  voltage  equal  to  the  rms  value  of  the  filtered 
signal.   The  rms  voltage  signal  was  recorded  as  a  function 
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of  time  on  an  X-Y  plotter.   The  system  output  was  a  graph 
of  the  scintillation  signal  instantaneous  relative  log- 
intensity  at  the  specific  frequency  as  a  function  of  time 
(Figure  8) .   The  graph  was  then  mechanically  integrated 
(using  a  planimeter)  to  find  the  average  relative  intensity 
of  the  signal  at  the  specific  frequency.   By  tuning  the 
bandpass  filter  to  a  variety  of  frequencies  and  recording 
the  signal  intensity,  the  scintillation  intensity  power 
spectrum  could  be  measured. 

The  tunable  bandpass  filter  had  a  low  frequency 
cut-off  at  20  Hz.   The  system  frequency  response  was  ex- 
tended to  5  Hz  by  playing  the  tape  of  the  scintillation 
signal  at  four  times  the  original  recording  speed  and  thereby 
increasing  the  original  signal  frequency  by  a  factor  of  four. 

A  tunable  oscillator  was  used  to  calibrate  the  wave 
analyzer  bandpass  filter  by  generating  a  pure  sinusoidal 
wave  at  the  desired  frequency  and  tuning  the  filter  to  match 
that  frequency.   An  oscilloscope  was  used  to  measure  the 
oscillator's  output  accurately. 

2.   Scintillation  Spectral  Analysis 

The  scintillation  power  spectrum  was  analyzed  to 
verify  that  the  system  did  not  distort  the  high  frequency 
components  of  the  scintillation  signal.   Spectrum  analysis 
of  three  one-minute  duration  signals  was  accomplished  using 
the  Fourier  spectrum  analysis  system.   The  three  signals 
analyzed  represented  the  maximum,  minimum  and  average  values 
of  C^Q,  recorded  during  the  first  Monterey  Bay  experiment. 
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The  relative  intensities  of  the  various  frequency  components 
were  averaged  to  give  a  representative  scintillation  spectrum 
as  shown  in  Graph  VIII. 

High  frequency  components  above  approximately  300  Hz 
had  negligible  power  levels.   The  system's  high  frequency 
distortion  limit  was  approximately  1500  Hz.   This  was 
determined  by  the  tape  recorder  frequency  response.   The 
experiment  showed  that  the  high  frequency  components  of  the 
scintillation  signal  were  not  distorted  by  the  system 
electronics . 

3.   Laser  Beam  Intensity  Profile  Analysis 

The  intensity  profile  of  the  laser  beam  was  examined 
to  determine  the  degree  of  truncation  by  the  output  optics. 
A  scanning  photomultiplier  was  used  to  measure  beam  intensity 
as  a  function  of  position  in  the  plane  of  the  output  aperture, 
Intensity  profiles  were  measured  with  and  without  the  output 
lens  in  position  and  examined  for  beam  truncation  by  the 
lens  holders. 

The  intensity  scans  showed  this  profile  to  be 
essentially  Gaussian  with  negligible  truncation.   The  pro- 
files also  indicated  that  the  laser  power  level  varied  in 
time,  during  the  warm-up  period.   Further  analysis  of 
intensity  stability  was  necessary  to  show  that  the  variations 
in  the  radiation  received  by  the  scintillometer  were  due  to 
atmospheric  propagation  effects  and  not  to  internal  intensity 
variations . 
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4.   Laser  Power  Stability  Analysis 

Laser  intensity  variations  were  monitored  and 
recorded  using  a  photomultiplier  tube  and  strip  chart 
recorder.   Measurements  of  intensity  as  a  function  of  time 
were  taken  in  the  center  and  wings  of  the  beam  profile. 
All  measurements  showed  large  fluctuations  in  beam  intensity 
during  the  warm-up  period.   These  fluctuations  were  in  the 
form  of  a  periodic  modulation  that  gradually  decreased  in 
amplitude  and  increased  in  period.   Maximum  modulation 
occurred  immediately  after  lasing  began.   After  approxi- 
mately one  hour  of  operation,  the  intensity  stabilized  and 
the  modulation  became  negligible.   However,  it  was  noticed 
that  a  5%  modulation  could  be  induced  in  the  beam  by  blowing 
air  across  the  laser  vents  after  the  laser  had  reached 
thermal  equilibrium.   The  cool  air  created  thermal  stresses 
on  the  laser  tube  and  mirror  supports,  and  caused  variations 
in  the  optical  gain. 

Realizing  that  during  the  second  Monterey  Bay  experi- 
ment, the  laser's  operating  environment  continuously  subjected 
the  laser  to  thermal  stress  from  meteorological  phenomena, 
the  validity  of  the  data  would  be  in  doubt  if  the  internal 
modulations  were  large.   To  test  the  effects  of  an  outdoor 
environment  on  the  intensity  stability,  the  laser,  photo- 
multiplier  and  strip  chart  recorder  were  assembled  on  the 
roof  of  Spanagel  Hall  and  the  experiment  was  repeated.   A 
short  propagation  path  minimized  atmospheric  scintillation 
in  the  beam. 
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The  laser  intensity  fluctuations  stabilized  from  35 
modulation  when  lasing  commenced,  to  approximately  4$ 
modulation  after  one  hour.   The  level  of  modulation  varied 
with  wind  velocity  but  never  exceeded  4%. 

Assuming  that  a  51  modulation  was  the  maximum 

internal  intensity  variation  encountered  during  the  second 

Monterey  Bay  experiment,  the  magnitude  of  the  perturbation 

effect  on  the  measured  value  of  C.^  was  calculated  from 

NO 


(59)  C^2  =  ([ln(I/(l})]2}   wh 


ere 


I   =   I   +  AI 
s     m 


In  the  equation,  I   is  the  instantaneous  intensity  of  the 

received  beam  due  to  scintillation,  and  AI  was  the  devia- 

*       m 

tion  from  the  average  intensity  due  to  internal  perturba- 
tions of  the  laser  output  intensity.   Calculations  using 
equation  (59)  with  values  of  AI   =  ±0.05(l)  and  a  typical 
probability  distribution  of  I   showed  that  the  maximum 

increase  in  C„0  values  was  6.4%  (computed  for  a  C^n  of 

-  8   -1/3 
2.0  x  10   m  '  ,  calculations  showed  that  the  perturbation 

effect  would  be  less  than  2.22%. 

These  calculations  represented  the  maximum  increase 

in  CNQ  values  due  to  a  continuous,  random  perturbation  of 

AI   =  ±0.05  (i)  .   In  actual  conditions,  the  measured  average 

intensity  perturbation  was  less  than  5%.   These  estimates 

are  an  upper  bound  for  possible  perturbation  effect  errors 
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in  the  data  taken  on  the  second  and  third  of  May.   The 
laser  was  shielded  inside  a  motel  room  during  the  first 
Monterey  Bay  experiment  and  perturbation  effects  should  have 
been  negligible.   Perturbation  effects  on  the  laser  used 
during  the  San  Nicolas  Island  experiment  by  the  research 
group  from  the  Naval  Missile  Center  are  an  unknown  quantity. 
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V.   DATA  ANALYSIS 

The  primary  goals  of  the  laser  propagation  research 
project  were: 

(1)  to  develop  instrumentation  and  expertise  for 
measurement  of  laser  beam  scintillation  and  micro- 
meteorological  turbulence, 

(2)  to  examine  the  correlation  between  micrometeoro- 
logical  atmospheric  turbulence  and  laser  beam 
scintillation  and 

(3)  to  verify  Tatarski's  [1]  formulation  of  the  theory 
of  laser  propagation  through  a  turbulent  atmosphere, 
in  the  marine  boundary  layer. 

Three  main  experiments  were  performed  between  January 
and  May  of  1973  during  which  optical  and  meteorological 
data  were  collected.   The  first  experiment  was  conducted 
on  January  10.   The  laser  radiation  was  transmitted  across 
Monterey  Bay  from  the  Holiday  Inn  to.  a  receiver  in  a  labora- 
tory at  Hopkins  Marine  Station  (see  Figure  9) .   Micrometeoro- 
logical  measurements  were  made  on-board  the  Naval  Postgraduate 
School's  Research  Vessel  Acania.   The  second  experiment  was 
conducted  at  San  Nicolas  Island  on  March  3.   A  laser  on  the 
Island  was  aimed  at  an  optical  detector  mounted  on  the  Acania. 
An  additional  set  of  laser  scintillation  measurements  was 
made  by  a  research  group  from  the  Naval  Missile  Center, 
Point  Mugu,  across  Laser  Bay  on  San  Nicolas  Island  (see 
Figure  10) .   The  third  experiment  was  conducted  across 
Monterey  Bay  between  the  Holiday  Inn  and  Hopkins  Marine 
Station  on  the  second  and  third  of  May.   Micrometeorological 
measurements  were  taken  on  the  beach  adjacent  to  the 
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Holiday  Inn.   The  data  collected  during  the  three  main 
experiments  will  be  reported  and  analyzed  in  this  section. 
General  conclusions  based  on  the  results  of  the  three 
experiments  will  be  presented  in  the  following  section. 

A.   FIRST  MONTEREY  BAY  EXPERIMENT 

The  first  Monterey  Bay  experiment  was  carried  out  on 
January  10.   The  laser  output  and  pointing  system  was 
placed  inside  a  room  in  the  Monterey  Holiday  Inn.   The  beam 
was  transmitted  a  distance  of  4.3  kilometers  across  the  Bay 
to  the  receiver  located  at  Hopkins  Marine  Station.   The 
transmission  path  was  horizontal  and  approximately  six  meters 
above  the  surface.   This  path  was  closely  transverse  to  the 
prevailing  wind  direction,  and  was  chosen  for  this  reason 
and  to  minimize  land  effects  on  turbulence  homogeneity.   The 
meteorological  data  were  taken  on-board  the  R.V.  Acania 
which  was  anchored  at  approximately  the  mid  point  of  the 
laser  transmission  path.   The  optical  and  micrometeorological 
measurements  were  taken  at  the  same  height  above  the  sea 
surface. 

The  optical  and  meteorological  estimates  of  CN  are 
listed  in  Data  Table  I  and  are  shown  as  functions  of  time  on 
Graphs  I  A  and  B.   Graph  II  is  a  scatter  plot  of  the  values 
of  CNQ  plotted  against  CNT« 

The  optical  measurements  of  C^  were  made  in  the 
laboratory  by  processing  the  tape  recorded  scintillation 
signal.   The  pulse  height  analyzer  accumulated  the  signal 
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probability  distribution  over  consecutive  five-minute 
periods.   This  corresponded  to  a  five-minute  integration  of 
the  laser  beam  scintillation  signal.   Meteorological  data 
were  recorded  on  magnetic  tape  and  later  transferred  to  a 
strip  chart  from  which  the  average  values  of  CT  were 
measured  for  five-minute  periods.   Lack  of  experience  in 
operating  the  scintillometer  system  caused  several  inter- 
ruptions in  data  collection. 

The  experimental  error  in  the  optical  measurements  of 
CNn  was  predominantly  due  to  errors  involved  in  measuring  the 

full  width  of  the  pulse  height  distribution.   For  values  of 

-  8  -1/3 
CNn  between  2  and  5  (X10   m  '  ),  the  standard  deviation  of 

the  measurements  was  approximately  5%.   Instrumentation 

error  and  the  error  involved  in  measuring  optical  wavelength 

and  optical  path  length  were  much  less  than  51. 

The  experimental  error  involved  in  the  meteorological 

measurements  of  CNT  was  predominantly  due  to  inaccuracies  in 

determining  the  average  value  of  CT  over  five-minute  intervals 

from  strip  chart  recordings  of  the  instantaneous  measurements 

of  CT.   To  complicate  this  process,  the  thermosonde  unit 

had  high  and  low  scales  for  different  extremes  of  turbulence. 

Whenever  the  scale  changed,  30  seconds  of  data  were  lost 

due  to  system  transients.   The  standard  deviation  was 

estimated  as  approximately  201  for  the  CNT  measurements. 

Errors  involved  in  measuring  pressure  and  temperature  and 

internal  system  errors  were  negligible  compared  to  averaging 

errors.   Experimental  error  bars  are  indicated  at  one  point 
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on  each  graph.   These  estimates  of  experimental  error 
apply  to  all  points  on  the  graph. 

In  analyzing  the  data,  it  is  important  to  realize  the 
fundamental  difference  between  the  nature  of  the  meteorologi- 
cal and  optical  measurements.   The  meteorological  measure- 
ments represent  the  degree  of  turbulence  at  only  one  point 
along  the  propagation  path.   The  optical  measurements 
represent  the  integration  of  turbulence  effects  on  laser 
propagation  over  the  total  propagation  path.   The  point 
measurements  of  turbulence  showed  greater  variations  in  the 
values  of  CNT  because  localized  cells  of  high  turbulence 
strongly  affected  the  meteorological  sensors.   The  optical 
measurements  did  not  respond  to  localized  cells  of  strong 
turbulence  unless  the  dimensions  of  those  cells  were  sig- 
nificant compared  to  the  propagation  path  length.   The  laser 
beam  scintillation,  in  effect,  averaged  the  effects  of 
atmospheric  turbulence  over  the  entire  optical  path.   This 
is  equivalent  to  saying  that  the  assumption  of  homogeneous 
turbulence  is  not  entirely  applicable  to  distances  on  the 
order  of  several  kilometers.   Hence,  One  should  analyze  the 
data  by  comparing  the  average  level  of  CNT  to  the  optical 

measurements  of  CXT~. 

NO 

The  data  showed  a  strong  correlation  between  the  optical 

and  meteorological  values  of  CN  with  the  exception  of  three 

data  regions  on  Graph  I.   At  approximately  1325,  the  value 

e  corre- 
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-  8         1/3 
of  CNT  increased  to  12.0  X  10   meters   '   but  the  corre 


sponding  value  of  CNQ  increased  to  only  4.36  x  10   meters 
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Again  at  approximately  1350,  the  value  of  CNT  increased 

- 8  -1/3 
to  15.0  x  10  m     whereas  CNn  remained  steady  at 

3.25  x  10~8m_1'3.   Finally,  between  1525  and  1630,  the 

values  of  CNT  showed  large  variations.   Corresponding  optical 

data  remained  relatively  constant.   The  average  value  of 

the  meteorological  data  was  approximately  equal  to  the 

optical  data. 

These  deviations  between  CMn  and  CNT  may  have  been 
caused  by  inhomogeneous  turbulence  which  strongly  affected 
the  temperature  sensors  but  was  not  sufficiently  wide- 
spread to  affect  the  laser  beam  propagation.   The  fact  that 
the  average  value  of  C„T  (between  1525  and  1630)  was  approxi- 
mately equal  to  the  relatively  steady  value  of  CN0  tended 
to  support  this  explanation.   However,  it  was  also  possible 
that  the  turbulence  was  homogeneous  and  that  the  optical 
values  of  CN  became  saturated.   Although  the  saturation 
effect  has  been  observed  in  over- land  measurements  by  Ochs 
in  1969  [3] ,  the  cause  of  saturation  has  not  been  clearly 
understood  or  explained. 

Graph  II  is  a  scatter  plot  of  the  optical  and 

meteorological  measurements  of  C„.      The  two  lines  drawn 

on  either  side  of  the  Cxrri-CXTrT,  equivalence  line  delineate  the 

NO   NT   n 

region  where  C^  equals  CNT  within  one  standard  deviation  of 
experimental  error.   Sixty-five  per  cent  of  the  data  points 
are  within  this  region  indicating  good  correlation  between 
measurements  of  CN0  and  CNT.   The  meteorological  values 
were  systematically  higher  than  the  optical  values.   This  ■ 
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would  be  expected  in  the  case  of  saturation.  However,  the 
lack  of  a  well-defined  saturation  level  for  optical  values 
precludes  a  definite  conclusion. 

This  experiment  indicated  that  the  optical  and 
meteorological  systems  that  had  been  developed  could 
measure  atmospheric  turbulence  and,  for  the  prevailing 
conditions  of  the  experiment,  the  local  thermal  fluctuations 
formed  a  good  predictor  of  beam  behavior. 

B.   SAN  NICOLAS  ISLAND  EXPERIMENT 

The  next  experiment  was  conducted  over  the  weekend  of 
March  3  at  San  Nicolas  Island.   Laser  beam  scintillation 
measurements  were  taken  over  a  one-kilometer  path  between 
San  Nicolas  Island  and  the  R.V.  Acania.   This  experiment 
tested  the  feasibility  of  making  scintillation  measurements 
between  the  shore  and  a  ship.   A  helium-neon  laser  on  the 
island  was  aimed  at  a  shipboard  detector  package.   Meteoro- 
logical measurements  were  made  on-board  the  Acania.   Inde- 
pendent measurements  of  CN0  were  made  by  the  optical 
propagation  group  from  the  Naval  Missile  Center,  Point  Mugu. 

The  data  from  this  experiment  are  listed  in  Data  Table 
II.   Graph  III  shows  the  measurement  o.f  CM  as  a  function  of 
time. 

The  problem  of  keeping  the  laser  beam  lined  up  with  the 
optical  receiver  proved  very  difficult.  The  laser  operator 
aimed  the  beam  at  a  retro-reflector  mounted  near  the  detec- 
tor and  used  the  reflection  for  sighting  corrections. 
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A  gimballed  mount  was  provided  for  the  detector  which  could 
be  aligned  with  the  laser  beam  and  thereby  reduce  the  effects 
of  ship  motion  on  the  detected  beam  intensity.   This  system 
proved  crude  and  awkward  to  operate.   By  making  measurements 
at  night,  the  need  for  a  background  compensating  photocell 
was  eliminated. 

The  major  problem  encountered  was  that  the  laser  beam 
could  not  be  pointed  at  the  detector  for  periods  of  time 
longer  than  five  to  ten  seconds  due  to  ship  motion.   The 
detector  was  moved  around  in  the  beam  by  the  roll  and  pitch 
of  the  ship  causing  a  kind  of  inverse  beam  wander  effect 
and  theoretically  increasing  the  measured  values  of  CNf.. 
Variations  in  laser  beam  intensity  were  due  to  both  scintil- 
lation and  ship's  motion.   However,  by  keeping  measurement 
times  very  short,  the  effect  of  ship  motion  was  minimized 
and  several  log-intensity  distributions  were  accumulated. 
Agreement  to  within  about  401  of  optical  measurements  along 
the  beach  line,  optical  measurements  from  shore  to  ship  and 
meteorological  measurements  on  shipboard  lead  to  some 
confidence  in  the  basic  measurement  techniques. 

This  experiment  showed,  however,  that  a  stabilized 
optical  system  was  necessary  to  keep  the  laser  beam  and 
detector  aligned.0  The  optical  system,  as  constructed,  could 
not  reduce  ship  motion  effects  to  an  acceptable  level.   It 
was  evident  that  the  system  needed  to  be  re-engineered  for 
shipboard  use. 
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C.   SECOND  MONTEREY  BAY  EXPERIMENT 

The  final  propagation  experiment  was  conducted  on  the 
second  and  third  of  May  across  Monterey  Bay.   As  before, 
the  laser  radiation  was  transmitted  from  the  Holiday  Inn 
to  Hopkins  Marine  Station  over  a  horizontal  path,  five 
meters  above  the  sea  surface.   The  optical  measurements 
were  taken  as  previously  described  for  the  first  Monterey  Bay 
experiment.   The  meteorological  measurements  were  taken  on 
the  beach  at  the  Holiday  Inn,  where  the  temperature  sensors 
were  positioned  approximately  two  meters  below  the  optical 
propagation  path.   The  values  of  CNT  were  extrapolated  for 
the  difference  in  measurement  levels. 

1.  Data  Taken  on  May  2,  19  73 

Administrative  and  technical  problems  caused  many 
interruptions  and  delays  in  data  accumulation  during  the  first 
day  of  the  experiment.   Numerous  beam  alignment  adjustments 
caused  the  system  to  be  shut  down  often  for  realignment. 
The  data  collected,  especially  before  1800,  were  of  ques- 
tionable accuracy.   Graphs  IV  A  and  B  show  the  optical  and 
meteorological  measurements  of  C^  plotted  as  functions  of 
time.   The  average  values  of  the  two  graphs  are  approximate- 
ly equal,  especially  for  data  taken  after  1800.   The  large 
variations  in  CNn  are  due  primarily  to  beam  alignment 
problems . 

2.  Data  Taken  on  May  3,  1973 

The  data  taken  on  the  second  day  of  the  experiment 
(Graphs  IV  C,  D,  E,  F,  G,  H)  showed  better  correlation 
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between  the  optical  and  meteorological  estimates  of  C„. 

Between  0830  and  1030,  the  values  of  CNT  increased  and 

varied  rapidly  in  time.   Corresponding  measurements  of 

C\t~  showed  a  lesser  increase  and  did  not  reach  the  high 
NO  & 

values  of  CNT.      This  was  probably  due  to  the  saturation 

of  CXTl~.   Local  cells  of  strong  turbulence  could  have  caused 
NO  6 

these  effects  but  it  was  doubtful  that  isolated  cells  would 

have  caused  sequential  increases  in  a  period  of  two  hours. 

Assuming  that  atmospheric  turbulence  over  the 

propagation  path  was  random,  one  expects  that  the  number  of 

localized  cells  of  atmospheric  turbulence  that  affect  one 

point  on  the  propagation  path  during  a  specific  time  period 

to  be  representative  of  the  number  of  cells  that  passed 

.any  point  along  the  path  during  the  same  time  period. 

There  was  no  reason  to  expect  that  localized  turbulence 

would  cluster  around  the  thermosonde  system.   One  or  two 

strong  increases  in  CNT  during  a  long  period  of  time  could 

be  blamed  on  isolated  localized  cells  of  high  turbulence 

that  would  have  little  effect  on  laser  beam  scintillation. 

However,  successive,  strong  increases  in  CNT  during  a 

long  period  of  time  could  only  be  due  to  an  overall 

increase  in  atmospheric  turbulence.   These  arguments  tend 

to  support  the  saturation  explanation  for  the  failure  of 

CXTrv   to   correlate  with  CXTrr 
NO  NT. 

From  1200  to  1600,  the  values  of  CNn  and  CNT 
showed  high  correlation,  while  atmospheric  turbulence 
remained  relatively  steady.   At  approximately  1630,  laser- 


59 


beam  alignment  was  lost  and  data  collection  ceased  until 
1845.   During  this  time  period,  the  values  of  C„T  showed  a 
gradual  decrease.   When  C^  measurements  were  re-commenced 
at  1845,  values  for  CN0  and  CNT  showed  a  high  degree  of 
correlation. 

3.   Data  Correlation 

Correlation  of  CNn  with  C„T  was  checked  by  plotting 
C^  as  a  function  of  CNT  in  Graph  V  (A,  B,  C)  and  Graph  VI 
(A,  B,  C,  D) .   Graph  V  A  shows  the  data  measured  on  the 
second  of  May.   Sixty-one  per  cent  of  all  the  data  points 
lie  within  one  standard  deviation  experimental  error  of  the 
C^NffC'NT  equivalence  line.   Experimental  error  bars  have  been 
plotted  on  only  one  point  but  apply  for  all  points.   No 
definite  level  of  saturation  was  observed.   Graphs  V,  B 
and  C,  show  the  data  from  the  second  of  May  sorted  according 
to  the  mean  wind  speed  at  the  time  of  measurement.   No 
definite  trends  were  noticed.   For  measurements  made  when 
the  mean  wind  was  less  than  seven  miles  per  hour,  66%  of 
all  the  data  points  lie  within  one  standard  deviation  of 
the  CNq-Cnt  equivalence  line.   For  measurements  made  when 
the  mean  wind  was  greater  than  seven  miles  per  hour,  59% 
of  all  data  points  lie  within  one  standard  deviation  of  the 
CNq-C„„  equivalence  line.   The  estimated  precision  and  the 
quality  of  the  data  and  the  lack  of  a  wide  range  of  wind 
speeds  do  not  justify  any  significant  conclusions  concerning 
the  effect  of  mean  wind  velocity  on  C^,  measurements. 
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Graph  VI  A  is  a  scatter  plot  of  the  C^  measurements 
made  on  the  third  of  May.   A  large  reduction  in  technical 
problems  led  to  a  higher  level  of  confidence  in  the  optical 
and  meteorological  data  taken  on  this  day.   Sixty-three  per 
cent  of  all  data  points  lie  within  one  standard  deviation 
of  the  CN0-CNT  equivalence  line.   Graphs  VI  B,  C,  and  D 
are  scatter  plots  of  CN  measurements  sorted  according  to  the 
mean  wind  velocity  at  the  time  of  measurement.   For  data 
taken  during  relatively  calm  periods  (Graph  VI  B) ,  data 
point  dispersion  is  very  large.   Only  231  of  all  the  points 
lie  within  one  standard  deviation  experimental  error  of  the 
equivalence  line.   The  data  taken  during  periods  when  the 
wind  speed  was  greater  than  three  but  less  than  six  miles 
per  hour  (Graph  VI,  C)  show  much  less  dispersion.   Seventy- 
nine  per  cent  of  the  data  points  are  located  within  one 
standard  deviation  experimental  error  of  the  CN0-CNT  equiva- 
lence line.   This  trend  toward  less  dispersion"  with  higher 
wind  speed  is  continued  in  Graph  VI   D  where  data  points 
measured  during  periods  when  the  wind  speed  was  greater 
than  six  miles  per  hour  are  plotted.   In  this  case,  931 
of  all  the  data  points  are  located  within  one  standard 
deviation  experimental  error  of  the  equivalence  line.   These 
data  indicate  that  an  increase  in  mean  wind  speed  was 
accompanied  by  an  increase  in  the  homogeneity  of  the  turbu- 
lence and  thus,  less  dispersion  of  data  points  on  the 
scatter  plots . 
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VI.   GENERAL  RESULTS 

Empirical  data  from  the  Monterey  Bay  experiment  indicate 

that  the  optical  and  micrometeorological  measurements  of 

atmospheric  turbulence  correlated  well  within  experimental 

tolerances.   Graph  VII  shows  that  681  of  all  measurements  of 

CXIO  and  CXTrr  were  within  one  standard  deviation  experimental 
NO      NT  r 

error  of  the  CXT„-CXT^  equivalence  line.   The  San  Nicolas 

NO  Nl   n 

Island  data  were  excluded  from  the  statistical  analysis 

because  of  the  unresolved  problem  of  ship  motion.   These 

results  indicate  that  it  is  possible  to  predict  laser  beam 

scintillation  at  a  wavelength  of  0.63  micrometers  in  the 

marine  boundary  layer  from  measurements  of  mean  pressure, 

mean  temperature,  and  the  covariance  of  the  temperature 

fluctuations . 

Some  evidence  indicates  that  saturation  of  CXT~ 

NO 

measurements  may  have  occurred  for  values  of  CN  greater 

-  8       - 1/3 
than  approximately  5.0  x  10   meters  '    .   However,  the 

lack  of  a  definite  level  of  saturation  in  the  scatter  plots 
suggests  that  definite  conclusions  on  this  must  await 
additional  measurements. 

The  correlation  between  the  optical  and  meteorological 
measurements  of  the  refractive  index  structure  constant 
increased  as  the  mean  wind  speed  increased.   Measurements 
taken  on  the  second  of  May  were  inconclusive  due  to  a  lack 
of  range  in  wind  speed  values.   It  was  expected  that   turbu- 
lence would  become  more  homogeneous  as  the  mean  wind  increased 
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The  theoretical  model,  from  which  the  relations  for  CN 
were  derived,  assumed  homogeneous  turbulence  in  the  hori- 
zontal plane.   Hence,  with  an  increase  in  wind  speed, 
turbulence  should  become  more  homogeneous  and  the  theoretical 
model  should  describe  the  physical  processes  more  accurately 
in  accordance  with  the  experimental  results. 

The  San  Nicolas  Island  experiment  demonstrated  the 
feasibility  of  making  scintillation  measurements  between 
the  R.  V.  Acania  and  a  shore  station.   However,  problems 
in  laser  alignment  due  to  ship  motion  indicated  that  the 
optical  receiver  system  needed  to  be  redesigned.   A 
stabilized  control  system  is  needed  to  keep  the  laser  beam 
and  detector  system  aligned. 

Finally,  the  high  correlation  between  the  optical  and 
meteorological  measurements  of  the  refractive  index  structure 
constant  indicates  that  Tatarski's  description  of  electro- 
magnetic propagation  is  essentially  valid  in  the  marine 
boundary  layer  at  a  wavelength  of  0.63  micrometers.   The 
validity  of  the  theory  appears  to  improve  with  increasing 
wind  speeds.   Additional  experiments  are  needed  to  obtain 
conclusive  evidence  to  support  Tatarski's  formulation  over 
a  wide  range  of  meteorological  conditions. 
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VII.   FUTURE  RESEARCH  SUGGESTIONS 

Refractive  index  structure  constant  measurements  should 
be  made  over  a  wide  range  of  turbulence  to  determine  the 
possible  existence  and  description  of  the  saturation  of 
optical  values.   These  measurements  could  further  substanti- 
ate or  refute  the  existing  results  of  this  project. 

Project  measurements  were  restricted  to  a  single 
wavelength  and  propagation  distance.   Additional  experiments 
should  be  made  to  determine  the  validity  of  the  theoretical 
model's  predicted  dependency  on  the  parameters  of  wavelength 
and  propagation  path  length.   Experiments  could  be  conducted 
at  various  wavelengths  using  various  types  of  lasers. 

Equipment  should  be  designed  to  shield  the  laser  from 
fluctuating  meteorological  conditions  and  to  continuously 
monitor  the  laser's  output  intensity.   This  would  reduce 
the  internal  intensity  modulation  and  provide  a  convenient 
method  of  determining  beam  stability. 

Additional  thermosonde  units  should  be  added  to  the 
meteorological  instrumentation  to  obtain  a  more  representa- 
tive sampling  of  turbulence  along  the  propagation  path. 
This  would  reduce  the  possibility  of  localized  isolated 
cells  of  inhomogeneous  turbulence  giving  false  indications 
of  total  path  turbulence.   Additional  experiments  should  be 
conducted  to  determine  if  the  shore  line  boundary  has  a 
significant  effect  on  the  meteorological  measurements. 
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Simultaneous  measurements  of  C,Trr  could  be  made  both  on  the 

NT 

beach  and  on-board  the  Acania  and  compared  for  correlation. 

A  digital  system  should  be  assembled  for  both  optical 
and  meteorological  data.   Digitizing  the  scintillation 
and  temperature  fluctuation  signals  would  greatly  reduce  the 
experimental  error  and  provide  a  direct  interface  with  the 
Naval  Postgraduate  School's  IBM  computer. 

Finally,  the  development  of  a  system  to  make  ship-to- 
shore  measurements  of  laser  beam  scintillation  would  be 
highly  desirable.   The  system  would  require  some  form  of 
stabilized  optics  to  reduce  ship  motion  effects.   This 
would  make  the  system  totally  portable  and  measurements 
could  be  made  almost  anywhere  on  the  coast. 
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VIII.  '  PROJECT  CONCLUSIONS 

The  experimental  measurements  obtained  during  this 
research  project  suggest  the  following  conclusions: 

(1)  Systems  were  developed  that  made  possible  simulta^ 
neous  optical  and  meteorological  measurements  of 
the  refractive  index  structure  constant  in  the 
marine  boundary  layer. 

(2)  Laser  beam  scintillation  at  0.63  micrometers  was 
predicted  from  meteorological  measurements  of 
atmospheric  temperature  fluctuations,  mean 
pressure  and  mean  temperature. 

(3)  Tatarski's  formulation  of  electromagnetic 
propagation  [1]  was  valid  for  propagation  in 
the  marine  boundary  layer  at  a  wavelength  of 
0.63  micrometers. 

(4)  The  validity  of  Tatarski's  formulation  increases 
with  increasing  mean  wind  speed  and  subsequent 
increasing  homogeneity  of  atmospheric  turbulence. 

(5)  More  research  is  needed  over  a  wider  range  of 
meteorological  conditions  to  substantiate  the 
generality  of  the  project's  conclusions.   These 
measurements  would  resolve  the  question  of  the 
possible  existence  and  description  of  a  satura- 
tion effect  on  optical  scintillation  measurements 
of  the  refractive  index  structure  constant. 
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Data  taken  on  the  second  of  May  over  Monterey  Bay. 


Sixty-one  per  cent  of  all  data  points  are  within  one 
standard  deviation  experimental  error  of  the  CN0-CNT 
equivalence  line. 
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LASER  OUTPUT  AND  POINTING  SYSTEM 
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Figure  2. 
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OPTICAL  DETECTOR  PACKAGE 
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Figure  3. 
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THE  LASER  SCINTILLOMETER 
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Figure  4 
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SCINTILLATION  SPECTRUM  ANALYZER 
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Figure  7. 
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APPENDIX  C 

DATA  TABLE  I 
FIRST  MONTEREY  BAY  EXPERIMENT 

Time  C^q  C^t 

1220  ---  3.21 

1225  2.87  2.80 

1230  2.50  3.02 

1235  ---  3.56 

1240  2.18  2.51 

1245  2.55  2.80 

1250  2.55  3.22 

1255  2.66  3.60 

1300  2.71  3.73 

1305  2.98  3.11 

1310  3.30  4.02 

1315  3.57  2.80 

1320  4.36  5.01 

1325  4.15  12.00 

1330  4.02 

1335  3.25  4.50 

1340  3.41  3.45 

1345  ---  3.53 

1350  3.25  15.00 

1355  ---  5.01 

1400  ---  3.02 

1405  3.00  2.03 

1410  2.87      '  2.52 

1415  2.93  3.73 

1420  2.77  3.54 

1425  ---  3.01 

1430  2.23   •  3.70 

1435  •  2.34       •  3.00 

1440  2.18  3.55 
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DATA  TABLE  I 
FIRST  MONTEREY  BAY  EXPERIMENT 

Time  C^-j.  C^T 

1445  2.18  3.11 

1450  2.87  3.52 

1455  2.55  2.91 

1500  2.55  2.43 

1505  2.55  2.72 

1510  2.39  2.84 

1515  2.18  2.45 

1520  •  2.39  2.61 

1525  2.55  4.10 

1530  2.55  3.04 

1535  2.55  1.40 

1540  2.61  2.61 

1545  2.61  3.51 

1550  ---  5.52 

1555  ---  4.51 

1600  2.66  1.54 

1605  2.71  3.40 

1610  .2.50  2.70 

1615  2.39  2.51 

1620  2.23  2.52 

1625         '  2.29  4.23 

1630  2.13  4.05 

1635  2.18 

1640  .  2.07 

1645  2.39 

1650  2.07 

1655  2.02 

1700  2.07 
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DATA  TABLE  III 

SECOND  MONTEREY  BAY  EXPERIMENT 

Time  (2  May)  CNQ  CNT 

1400  ---  15.0 

1405  ---  14.6 

1410  6.01  14.6 

1415  ---  15.0 

1420  5.60  14.6 

1425  ---  2.67 

1430  5.52  3.46 

1435  5.60  3.46 

1440  5.79  3.58 

1445  5.60  2.13 
1450 

1455  5.01 

1500  5.34 

1505  5.25 

1510  6.21 
1515 

1520  ---  "  --- 

1525  ---  5.94 

1530  ---  4.81 

1535  ---  6.13 

1540  ---  5.33 

1545  ---  5.12 

1550  5.01  2.67 

1555  5.69  6.94 

1600  5.90  5.26 
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DATA  TABLE  III 
SECOND  MONTEREY  BAY  EXPERIMENT 


Time  (2  May) 

CNO 

1605 

5.31 

1610 

7.01 

1615 



1620 

6.10 

1625 

6.72 

1630 

5.71 

1635 



1640 

6.25 

1645 



1650 

6.39 

1655 



1700 



1705 



1710 



1715 



1720 



1725 

5.70 

1730 



1735 



1740 



1745 

4.91 

1750 

4.07 

1755 

4.35 

1800 

4.92 

CNT 


5.33 
5.12 
4.61 
4.91 
4.29 
4.61 
4.61 
5.12 
4.30 
4.91 
4.81 
4.71 
5.12 
4.40 
5.02 
4.91 
4.81 
4.61 
4.74 
4.40 
4.71 
4.81 
4.30 
4.39 
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DATA  TABLE  III 
SECOND  MONTEREY  BAY  EXPERIMENT 


Time  (3  May)        CNQ  CNT 


1805  5.35  4.59 

1810  ---  4.59 

1815  4.47  4.69 

1820  4.77  4.59 

1825  4.07  4.69 

1830  4.49  4.69 

1835  3.79  4.39 

1840  3.51  4.59 

1845  3.23  4.49 

1850  3.65  4.39 

1855  3.79  4.59 

1900  4.57  4.69 

1905  4.07  4.48 

1910  4.07  4.28 

1915  4.07  4.51 

1920  3.93  4.48 

1925  4.77  4.18 

1930  4.50 
1935 

1940  3.79 

1945  3.65 

1950  3.65 
1955 
2000 
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DATA  TABLE  III 
SECOND  MONTEREY  BAY  EXPERIMENT 

Time  (3  May)  CNQ               CNT 

0630  2.72 

0635 

0640  2.25 

0645 

0650 

0655  2.53 

0700  2.53 

0705  2.34 

0710  2.06 

0715  2.06 

0720  1.59 

0725  1.97 

0730  2.34               2.31 

0735  2.34               3.08 

0740  2.15               1.95 

0745  2.15             .  1.64 

0750  2.25               3.60 

0755  2.62               3.60 

0800  2.62               4.01 

0805  2.28               3.70 

0810  2.15               3.17 

0815  2.62               3.70 

0820  2.81               3.91 

0825  2.72               4.63 

0830  2.72               5.14 

0835  3.47               6.16 

0840  3.56               9.37 

0845  3.56               5.45 

0850  3.56               8.03 
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DATA  TABLE  III 
SECOND  MONTEREY  BAY  EXPERIMENT 

Time  (3  May)  CNQ  CNT 

0855  3.84  6.49 

0900  4.31  5.81 

0905  4.12  5.10 

0910  4.50  6.12 

0915  .  4.31  4.90 

0920  4.22  7.97 

0925  3.84  6.12 

0930  3.84  6.12 

0935  ---  5.10 

0940  ---  4.39 

0945  ---  4.49 

0950  2.35  5.10 

0955  2.97  7.97 

1000  2.15  5.61 

1005  ---  3.82 

1010  3.09  4.25 

1015  3.09  4.25 

1020  3.28  8.85 

1025  3.28  8.85 

1030  3.47  4.70 

1035  ---  4.27 

1040  ---  3.50 

1045  4.22     •  4.19 

1050  4.22  3.84 

1055  3.71  3.93 

1100  3.56  3.42 
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DATA  TABLE  III 
SECOND  MONTEREY  BAY  EXPERIMENT 


Time  (3  May) 


CNO 

CNT 

3.47 

3.50 

3.47 

3.67 

3.47 

3.75 

5.15 

3.84 

5.90 

3.07 



3.59 

4.12 

3.95 

4.03 

3.85 

4.12 

3.85 

4.87 

4.14 

4.78 

4.14 

4.31 

3.73 

4.22 

4.11 

3.94 

3.92 

4.40 

4.20 

4.03 

4.11 

4.69 

4.39 

4.69 

4.20 

5.45 

3.83 

5.70 

4.48 

5.55 

4.67 

4.31 

4.67 

4.50 

4.47 

4.03 

4.29 

1105 
1110 
1115 
1120 
1125 
1130 
1135 
1140 
1145 
1150 
1155 
1200 
1205 
1210 
1215 
1220 
1225 
1230 
1235 
1240 
1245 
1250 
1255 
1300 


99 


DATA  TABLE  III 
SECOND  MONTEREY  BAY  EXPERIMENT 


Time  (3  May) 

.   CNO 

1305 

4.03 

1310 

4.03 

1315 

4.12 

1320 

4.31 

1325 

4.40 

1330 

4.59 

1335 

4.50 

1340 

4.78 

1345 

4.12 

1350 

3.65 

1355 

3.47 

1400 

3.47 

1405 

3.56 

1410 

3.56 

1415 

3.56 

1420 

3.56 

1425 

3.92 

1430 

3.75 

1435 

3.47 

1440 

3.56 

1445 



1450 



1455 

3.34 

1500 

3.75 

r 

NT 


3.92 
3.73 
4.20 
4.29 
4.39 
4.50 
4.69 
4.41 
4.41 
4.69 
4.31 
4.50 
4.50 
4.50 
4.41 
4.78 
4.69 
4.50 
4.50 
4.60 
4.50 
4.22 
3.94 
4.22 
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DATA  TABLE  III 
SECOND  MONTEREY  BAY  EXPERIMENT 


Time  (3  May) 

LNO 

1505 

4.12 

1510 

4.03 

1515 

4.03 

1520 

4.03 

1525 

5.25 

1530 



1535 



1540 

4.40 

1545 



1550 



1555 

5.15 

1600 

5.05 

1605 

4.69 

1610 

3.75 

1615 

4.40 

1620 



1625 

3.75 

1630 

4.12 

1635 



1640 



1645 



1650 



1655 



1700 



CNT 


4.31 
3.94 
4.12 
3.75 
3.94 
4.41 
4.22 
4.13 
4.22 
4.22 
3.94 
4.12 
3.94 
4.31 
3.94 
3.75 
3.94 
4.24 
4.52 
2.82 
2.73 
3.18 
2.63 
2.82 
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DATA  TABLE  III 
SECOND  MONTEREY  BAY  EXPERIMENT 

Time  (3  May)  CNQ  CNT 

1705  ---  3.01 

1710  ---  2.35 

1715  ---  2.07 

1720  ---  2.26 

1725  ---  2.35 

1730  ---  2.48 

1735  ---  2.48 

1740  ---  2.48 

1745  ---  1.88 

1750  ---  1.88 

1755  ---  2.06 

1800  ---  2.15 

1805  ---  2.24 

1810  ---  2.05 

1815  .  ---  1.86 

1820  ---  1.96 

1825  ---  1.86 

1830  ---  1.96 

1835     '  ---  1.68 

1840  1.97  1.12 

1845  1.87  1.68 

1850  1.78  1.31 

1855  1.68             •  2.01 

1900  1.87  2.05 


102 


DATA  TABLE  III 
SECOND  MONTEREY  BAY  EXPERIMENT 

Time  (3  May)  CNQ  CNT 

1905  1.97  2.34 

1910  2.06  2.34 

1915  2.06  1.96 

1920  2.25            •   2.38 

1925  2.53  2.01 

1930  2.72  1.87 

1935  2.44  1.68 

1940  2.53  1.22 

1945  2.06  1.40 

1950  2.25  1.78 

1955  2.34  1.80 

2000  2.34  1.68 

2005  1.97  1.73 

2010  1.78  1.31 

2015  1.78  1.78 

2020  1.68  1.63 

2025  1.87  1.78 

2030     .  2.06  1.63 

2035  1.97  1.78 

2040  1.87  1.87 

2045  2.06  2.10 

2050  1.96  2.29 

2055  2.06 

2100  1.69 
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